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Abstract 

We construct new semi-realistic Type IIB flux vacua on Z2 X Z2 orientifolds with three- and 
four- Standard Model (SM) families and up to three units of quantized flux. The open-string 
sector is comprised of magnetized D-branes and is T-dual to supersymmetric intersecting D6- 
brane constructions. The SM sector contains magnetized D9-branes with negative D3-brane charge 
contribution. There are large classes of such models and we present explicit constructions for 
representative ones. In addition to models with one and two units of quantized flux, we also 
construct the first three- and four-family Standard-like models with supersymmetric fluxes, i.e. 
comprising three units of quantized flux. Supergravity fluxes are due to the self-dual NS-NS and 
R-R three- form field strength and they fix the toroidal complex structure moduli and the dilaton. 
The supersymmetry conditions for the D-brane sector fix in some models all three toroidal Kahler 
moduli. We also provide examples where toroidal Kahler moduli are fixed by strong gauge dynamics 
on the "hidden sector" D7-brane. Most of the models possess Higgs doublet pairs with Yukawa 
couplings that can generate masses for quarks and leptons. The models have (mainly right-) chiral 
exotics. 



I. INTRODUCTION 



One of the challenging and essential problems in string theory is the construction of 
realistic string vacua, which can stabilize the moduli fields, generate Standard Model (SM)- 
like gauge structure and induce a (de Sitter) cosmological constant with supersymmetry 
breaking. Such constructions would provide a bridge between string theory and realistic 
particle physics. M-theory provides a framework where, in addition to perturbative heterotic 
string vacua, the physical string vacua could be probed in the perturbative Type I, Type 
IIA and Type IIB superstring theory. In particular, the discovery of D-brane dynamics 
makes it possible to construct consistent four-dimensional super symmetric N — 1 chiral 
models with non-Abelian gauge symmetry on Type II orientifolds, by employing conformal 
field theory techniques in the open string sector. The first such supersymmetric models 
were based on Z2 x Z2 orientifolds [1, 2]. [Non-supersymmetric constructions were given 
in [3, 4, 5, 6] (see also [7] and for earlier work [8, 9]).] Subsequently, a number of SM- 
like models, GUT models and their variations have been constructed in various orbifold 
backgrounds, and the associated phenomenology has been discussed. [For a partial list of 
non-supersymmetric constructions, see [10], further supersymmetric constructions are given 
in [11]- [14] and further developments in connection with the study of effective couplings and 
phenomenological implications, see [15] -[21] and references therein.] 

[Recently important progress has been made in constructions of supersymmetric chiral 
solutions of Type II Gepner models; see [22, 23] and references therein. Specifically, the 
recent impressive results of [23] provide large classes of three-family Standard-like Models 
with no chiral exotics. Note however, that these exact conformal field theory models are 
located at the special points in the moduli space where the geometric picture is lost. In 
particular couplings, such as Yukawa couplings, do not possess hierarchies associated with 
the size of the internal spaces, such as in the case of the toroidal orbifolds with D-branes. In 
addition, due to the lack of geometric interpretation, the introduction of supergravity fluxes 
does not seem to be possible. ] 

In spite of these successes, the moduli stabilization in open string and closed string sectors 
remained an open problem, even though in some cases some complex structure moduli (in 
the Type IIA picture) and dilaton fields may be stabilized due to non-perturbative gauge 
dynamics, associated with the gaugino condensation in the hidden sector (see, e.g., [24].) 
Turning on supergravity fluxes introduces a supergravity potential, which provides another 
way to stabilize the compactification moduli fields by lifting continuous moduli space of 
the string vacua in the effective four- dimensional theory (see, e.g., [25]). However, the 
introduction of supergravity fluxes imposes strong constraints on consistent constructions, 
since such fluxes modify the global Ramond-Ramond (RR) tadpole cancellation conditions. 
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Meanwhile, the fluxes will typically generate a back-reaction on the original geometry of the 
internal space, thus changing the nature of the internal space. 

On the Type IIA side the supersymmetry conditions of flux compactifications are less 
understood. Nevertheless recent work [26, 27] revealed the existence of unique flux vacua 
for massive Type IIA string theory with SU(3) structure, whose geometry of the internal 
six-dimensional space is nearly-Kahler and four- dimensional space is anti-de Sitter (AdS) 
(for the discussion on necessary and sufficient conditions of N=l compactifications of mas- 
sive IIA supergravity to AdS (4) with SU(3) structure, see also [28]). One such example is 
the S ^2) — ^3 x ^3 cose t space that has three supersymmetric three-cycles that add up 
to zero in homology [27, 29]. Therefore the total charge of the D6-branes wrapping such 
cycles is zero and no introduction of orientifold planes on such spaces is needed. Moreover, 
since the three-cycles intersect pair-wise, the massless chiral matters appear at these inter- 
sections. This construction [27] therefore provides an explicit example of supersymmetric 
flux compactifications with intersecting D6-branes. Further progress has also been made 
in the construction of N = 1 supersymmetric Type IIA flux vacua with SU(2) structures 
[27], leading to examples with the internal space conformally Calabi-Yau. However, explicit 
constructions of models with intersecting probe D6-branes for such flux compactifications is 
still awaiting further study. 

On the Type IIB side the intersecting D6-brane constructions correspond to models with 
magnetized branes with the role of the intersecting angles played by the magnetic fluxes on 
the branes. The dictionary for the consistency and supersymmetry conditions between the 
two T-dual constructions is straightforward, see e.g., [30, 31]. The supersymmetric Type IIB 
flux compactifications are also better understood; see, e.g., [32]-[35], [30, 31] and references 
therein. In particular, examples of supersymmetric fluxes and the internal space conformally 
Calabi-Yau are well known. The prototype example is a self-dual combination of the Neveu- 
Schwarz-Neveu-Schwarz (NS-NS) H z and RR F 3 three-forms, corresponding to the primitive 
(2,1) form on Calabi-Yau space. Since the back- reaction of such flux configurations is mild, 
i.e., the internal space remains conformal to Calabi-Yau, these Type IIB flux compactifica- 
tions are especially suitable for adding the probe magnetized D-branes in this background. 
However, the quantization conditions on fluxes and the modified tadpole cancellation condi- 
tions constrain the possible D-brane configurations severely. In Refs. [30, 31] the techniques 
for consistent chiral flux compactifications on orbifolds were developed, however, no explicit 
supersymmetric chiral SM-like models were obtained. 

Most recently, by introducing magnetized D9-branes carrying negative D3-brane charges 
in the hidden sector, in Ref. [36] the first example of three-family SM-like string vacuum with 
one unit of quantized flux turned on was obtained, and subsequently, the first four-family SM- 
like string vacuum with one unit of fluxes was constructed in [37]. These constructions could 
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be T-dual to the supersymmetric models of intersecting D6-branes on Z 2 x Z 2 orientifold 
with the Sp(2) L x Sp(2) R or Sp(2f) L x Sp(2f) R gauge symmetry in the electroweak sector, 
respectively [14, 18]. [Without fluxes, the first models of that type were toroidal models 
with intersecting D6-branes [18] where the RR tadpoles were not explicitly cancelled. (For 
the subsequent generalization to tilted tori see [38].) The Z 2 x Z 2 orientifold construction 
in [14] provided the first model of that type that cancelled RR-tadpoles by introducing 
an additional stack of D6-branes with unitary symmetry; in the T-dual picture those are 
precisely the magnetized D9-branes with the negative contribution to the D3-brane charge.] 
In spite of these successes, we are confronted by a number of problems: 

(i) The semi-realistic SM-like string vacua with four-dimensional N — 1 supersym- 
metric fluxes have not been constructed, yet. In view of this drawback, effects of non- 
super symmetric fluxes, as a key mechanism for breaking supersymmetry, has been addressed 
[39, 40, 41]. This analysis [39] leads to soft supersymmetry breaking masses M so f t ~ 
where M s and M P i are the string scale and Planck scale, respectively. This result implies 
an intermediate string scale or one has to introduce an inhomogeneous warp factor in the 
internal space in order to stabilize the electroweak scale. 

(ii) These flux vacua stabilize the dilaton and toroidal complex structure moduli. How- 
ever, the Kahler moduli do not enter the flux-induced superpotential and hence are hard to 
be completely fixed. So, we are still typically faced with the vacuum degeneracy problem. 
[The Kahler moduli fields in Type IIB string theory are T-dual to the complex structure 
moduli in Type IIA string theory (intersecting D6-brane scenario). In the T-dual picture, 
the latter moduli can often be stabilized by employing non-perturbative dynamical mech- 
anism, such as the gaugino and matter condensation in the hidden sector. However, these 
mechanisms are difficult to employ on the Type IIB side due to the additional matter content 
on the associated magnetized D-branes.] 

(iii) For explicit Type IIB orientifolds the imaginary self-dual fluxes are quantized in 
rather large flux units, e.g., for Z 2 x Z 2 the orientifolds elementary flux unit is 64. Therefore 
the constructions of semi-realistic flux vacua is very constrained; only one unit of flux is al- 
lowed for known semi-realistic three- [36] and four-family [37] models. Thus the introduction 
of fluxes is restricted to very few known semi-realistic examples and is not typical. 

Following our previous work [37], in this paper we systematically study new constructions 
of three- and four-family SM-like string vacua with supergravity fluxes on Type IIB Z 2 x Z 2 
orientifolds. The major technical difficulty in constructions of semi-realistic flux vacua on 
Type IIB orientifolds is to ensure the cancellation of the large positive D3-brane charge 
contribution to RR tadpoles by the fluxes. Similar to the D-brane models without fluxes 
[14], and the subsequent work with fluxes [36, 37], the important role in tadpole cancellation 
is played (in the Type IIB picture) by magnetized D9-branes which carry large negative 
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D3-brane charges. In the past constructions ([14, 36, 37]) such D9-branes were introduced 
as a part of the "hidden sector" . 

In this paper, we consider new types of constructions where the magnetized D9-branes 
with large negative D3-brane charges are introduced as a part of the SM sector. For this new 
setup, we find that the constructions of SM-like flux vacua are much less constrained and 
obtained a large class of new models. In particular, in addition to many new models with 
one unit of quantized flux, we obtain first three- and four-family models with two units of 
quantized charge, as well as the first three- and four-family examples with super symmetric 
flux, i.e. three units of quantized flux. Such supersymmetric three- and four-family SM-like 
models have toroidal Kahler moduli fixed by supersymmetry conditions and the string scale 
can be close to the Planck scale. These models have (mainly right-) chiral exotics. However, 
most of the models have Higgs doublet pairs with Yukawa couplings to quarks and leptons 
and thus may generate the SM fermion mass hierarchies, explain the CKM quark mixing 
matrix and the PMNS neutrino mixing matrix at the tree level, and even give large masses 
to some of the bi-fundamental chiral exotics. Finally, with the open-string moduli fixed by 
the flux-induced soft masses, we are able to construct the first SM-like string vacua with 
strong infrared gauge dynamics on the "hidden sector" D7-branes, which leads to gaugino 
condensation and provides examples where the third toroidal Kahler modulus is fixed by 
the strong gauge dynamics a la KKLT [45]. 

The paper is organized in the following way. In Section II, we systematize the con- 
structions of supersymmetric string vacua with supergravity fluxes on Type IIB Z 2 x Z 2 
orientifolds. In Section III, we classify the classes of the SM-like flux vacua on Type IIB 
Z 2 x Z 2 orientifolds. Subsequently, we discuss in detail explicit constructions of the repre- 
sentative models with one, two and three (supersymmetric) units of quantized flux, as well 
as a specific construction with gaugino condensation on the "hidden sector" D7-branes. In 
the Appendix we provide tables of all the explicitly constructed representative models. We 
conclude with discussions and open problems in Section IV. 

II. MAGNETIZED D-BRANES AND TYPE IIB FLUX COMPACTIFICATIONS 

ON T 6 /(Z 2 x Z 2 ) ORIENTIFOLDS 

Flux compactifications on simplest toroidal T 6 Type IIB orientifolds, on which most of 
the previous work has focused, are unlikely to provide a framework for constructions of 
semi-realistic flux vacua. We shall therefore focus on the simplest orbifold constructions, 
i.e. on T 6 /(Z 2 x Z 2 ) orientifolds. 

The internal space T 6 is chosen to be factorized as a direct product of three two-tori, i.e. 
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T 6 = T 2 x T 2 x T 2 , whose complex coordinates are z i: % — 1, 2, 3 for the i-th two-torus, 
respectively. The generators 9 and uo for the orbifold group Z 2 x Z 2 , act on the complex 
coordinates of T 6 as 



9: (z 1 ,z 2 ,z 3 ) ^ (-z 1 ,-z 2 ,z 3 ) , 
UO : (Z!, z 2 , z 3 ) -> (zi, -z 2 , -2 3 ) • 



(1) 



On Z 2 x Z 2 orbifold Type IIB string theory contains in the untwisted sector the four- 
dimensional N = 2 supergravity multiplet, the dilaton hypermultiplet, hu hypermultiplets 
and h 2 i vector multiplets, which are all massless. The orbifold action projects out several 
components of the metric of a general T 6 geometry and, as a result, we are left with fewer 
Kahler and complex structure parameters. These are encoded for the untwisted moduli in 
terms of the Hodge numbers, as (h n , h 21 ) unt = (3, 3). On the other hand, each of the three 
elements 9, uo and 9uo has a fixed-point set given by 16 T 2 's, and the corresponding twisted 
sectors also contribute to the Hodge numbers of the orbifold. For a particular choice of 
discrete torsion, this contribution is given by (hu, h 2 i) tw = (0, 3 x 16). The contributions 
from both, the untwisted and twisted sectors hence add up to (hu, h 21 ) = (3, 51). 

Orientifold planes are necessary for the introduction of the open-string sector, and the 
associated orientifold projection can be denoted by QR, where Q is the world-sheet parity 
projection and R (acting on Type IIA as the holomorphic Z 2 involution) acts on the complex 
coordinates as: 



Thus, the model contains 64 03-planes and 4 07;-planes, transverse to the % two-torus. 
This orientifold action projects the above N = 2 spectrum to an iV = 1 supergravity 
multiplet, the dilaton chiral multiplet, and 6 untwisted and 48 twisted geometrical chiral 
multiplets. 

In order to cancel the negative RR charge contributions, due to these O-planes, we need to 
introduce D(3+2n)-branes which are filling up four-dimensional Minkowski space-time and 
wrapping 2n-cycles on the compact manifold. We choose the construction with magnetized 
D-branes. [A detailed discussion for toroidal/orbifold compactifications with magnetized 
D-branes is given, e.g., in [30].] Concretely, for one stack of N a D-branes wrapped m % a times 
on the i th two-torus T 2 , we turn on n l a units of magnetic fluxes F a for the center of mass 
U(l) a gauge factor on each T 2 , such that 



Hence, the topological information of this stack of D-branes is encoded in iV a -number of D- 
branes and the co-prime number pairs (n l a , m % a ). The D9-, D7-, D5- and D3-branes contain 0, 



R : (zi, z 2 , z 3 ) -> (-z ± , -z 2 , -z 3 ) . 



(2) 




(3) 
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1, 2 and 3 vanishing m l a s, respectively. Introducing for the i th two-torus the even homology 
classes [Oj] and [T 2 ] for the point and the two-torus, respectively, the vectors of RR charges 
of the a th stack of D-branes and its image are 

[n a ] = n km + ™M)> = nli «m - ™m) , (4) 

i=i 

respectively. Similarly, for the 03- and 07j-planes appearing on T 6 /(Z 2 x Z 2 ) orientifold 
which respectively correspond to ilR, VlRuo, VLROuo and VLR6 O-planes, we have 

tlR : [n 03 ] = [Oi] x [0 2 ] x [0 3 ]| 

QRj: [n D7l ] = -[Oi] x [7f] x [T 2 ]; 

fii^w: [n D72 ] = -[T 2 ] x [0 2 ] x [T 3 2 ]; 

fii?^: [n O73 ] = -[T 1 2 ]x[T 2 2 ]x[0 3 ]. (5) 

The "intersection numbers", which determine the chiral massless spectrum, are 

3 

I ab = [U a ] ■ [U b ] = HKml - v? b m\) , 
i=i 

3 

I aV = [ILJ • [Yi v ] = - JJ«mj + njm* ) , 

i=i 
3 

/ aa , = [n a ]-[n a; ] = -8n«0 , 

i=l 

J a0 = ^[n a ] • [II 0p ] = 8(-m> 2 m3 + m i n 2 n 3 + n i m 2 n 3 + n i n 2 m 3) ^ (6) 
p 

where [n 0p ] = [n 03 ] + [n 07l ] + [n 072 ] + [Il 073 ] is the sum of 03-plane and 07 r plane homology 
classes. Similar to the discussions in [2], the physical chiral spectrum should be invariant 
under the full orientifold symmetry group and is tabulated in Table I. Flux vacua on Type 
IIB orientifolds with four-dimensional N — 1 supersymmetry are primarily constrained by 
the RR tadpole cancellation conditions and conditions for iV = 1 supersymmetry in four 
dimension, which we describe in detail in the following subsections. 

A. RR Tadpole Cancellation Conditions 

In the Type IIB picture the fluxes, we consider, are due to the self-dual three-from field 
strength, which contributes to the D3-brane field strength equation of motion, and thus 
modifies the D3 charge conservation (RR-tadpole cancellation) conditions on the compact 



7 



TABLE I: General spectrum for magnetized D-branes on the Type IIB T 6 jiZi x Z 2 ) orientifold. 
The representations in the table refer to U(N a /2), the resulting gauge symmetry due to Z 2 x Z 2 
orbifold projection. For supersymmetric constructions, scalars combine with fermions to form 
chiral supermultiplets. 



Sector 


Representation 


aa 


U(N a /2) vector multiplet 
3 adjoint chiral multiplets 


ab + ba 


lab (Pa,Ob) fermions 


ab' + b'a 


lab' (Ha, Ob) fermions 


aa' + a'a 


\ (Jaa' ~ \la,Op) HJ fermions 
\ (ha' + \la,o P ) B fermions 



orientifold. The RR charges, carried by magnetized D-branes, are classified by their asso- 
ciated homology classes. Explicitly, for one stack of N a D-branes with wrapping numbers 
(774,7774), it carries D3-, D5-, D7- and D9-brane RR charges 

Q3 a = N a nln 2 a n 3 a , (Q5 t ) a = N a rrt a n{n k a , 
(Q7i) a = N a nimim k a , Q9 a = N a m l a m 2 a ml, (7) 

where % 7^ j 7^ k, and a permutation is implied for (Q5i) a and (Q7i) a - So, the RR tadpole 
cancellation conditions can be described as 

£ N a [n o ] + N a t n a'] + £ N 0p Q p l^Op] + ^flux = , (8) 

a a p 

where the third term contribution comes from the 03- and 07j-planes, with No p and Qo P 
denoting their numbers and RR charges, respectively. And A^a ux is the amount of the fluxes 
turned on, and is quantized in units of the elementary flux. 

For a supersymmetric Dp/Dp'-brane system on Type IIB T 6 /(Z 2 x Z 2 ) orientifold, only 
D3- and D7-branes are allowed to be wrapped along the orientifold planes. [In the following 
we shall refer to this type of branes as "filler branes" .] Given that N 0p Q 0p = 2 9 ~ p (— 2 P ~ 4 ) = 
—32 in Dp-brane units for Sp-type O-planes, the RR tadpole cancellation conditions can be 
further simplified as 

_Ar(o)_^g 3a _ ijV flux = -16, 
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- N ® + J2(Q7i)a = -16, i*j*k, (9) 

a 

where and with 2 =1, 2 and 3 respectively denote the number of filler branes, i.e. D- 
branes which wrap along the 03- and 07j-planes and only contribute to one of the four kinds 
of D3- and D7-brane charges. As for D5- and D9-brane RR tadpoles, their cancellations are 
automatic since these D-branes and their QR images carry the same absolute value of the 
corresponding charges, but with opposite sign of charges. 



B. Conditions for Four-Dimensional N = 1 Supersymmetry 

Four-dimensional N — 1 supersymmetric vacua from flux compactification require that 
1/4 supercharges of the ten-dimensional (T-dual) Type I theory be preserved in both open 
string and close string sectors. We shall discuss both sectors separately. 

For the closed string sector, the specific Type IIB flux solution on orientifolds comprises 
of self-dual three- form field strength and it has been discussed, e.g., in [32, 44]. While RR 
F 3 and NSNS H 3 three-form fluxes are turned on, the induced three-form G 3 = F 3 — tH 3 , 
with r = a + i/g s being Type IIB axion-dilaton coupling, contributes to the D3-brane RR 
charges 

AW = (A ] ,, 2 / H 3 AF 3 = I -L / G3AG3, (10) 
(4vrV) 2 J X(j (4ttV) 2 2tj J Xf> 

where 77 is the imaginary part of the complex coupling r. Dirac quantization conditions 
of F 3 and H 3 on T 6 /(Z 2 x Z 2 ) orientifold require that iV flux be a multiple of 64, and the 
BPS-like self-duality condition: * 6 G3 = iG 3 ensures that its contribution to the RR charges 
is positive. Supersymmetric configuration implies that G 3 background field should be a 
primitive self-dual (2,1) form. A specific supersymmetric solution which is useful for our 
purpose is [30, 44] 

G 3 = -j= e -7 ™/ 6 ( dzidz 2 dz 3 + dz\dz 2 dz 3 + dz\dz 2 dz 3 ), (11) 

where the additional factor 4 is due to the Z 2 x Z 2 orbifold symmetry. The fluxes stabilize 
the complex structure toroidal moduli at values 

ri = r 2 = r 3 = r = e 2 ™/ 3 , (12) 

leading to the RR tadpole contribution in Eq. (9) 

AW = 192. (13) 
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This result therefore implies that in order to construct supersymmetric SM-like flux vacua, 
we have to introduce at least the D3 charge conservation condition, which is thus hard to 
achieve. 

In the open-string sector, for D-branes with world- volume magnetic field F % = 
the four-dimensional N — 1 supersymmetry is ensured if and only if ^ 6i = (mod 2tt) 
is satisfied [30]. Here the " angle" 0j with the range {0, 2tc}, is determined in terms of 
the world-volume magnetic field as tan(^) = (F 1 ) -1 = I ^ and \ % — R\R\i the area 
of the % th two-torus Tf in a' units, is the Kahler modulus of the i th two-torus Tf. This 
supersymmetry condition can be cast in the form: Yliii^ 1 ) 1 ~ (F 1 F 2 F 3 )^ 1 = 0, along 
with Ei<j {F'Fiy 1 - 1 < for ni n,n k > or ^ t<J (F l F j ) _1 - 1 > for ni n,n k < 0, 
which can be rewritten in the following form: 

- x A Q9 a + x B (Q5i) a + x c (Q5 2 ) a + x D (Q5 3 ) a = 0, 

-Q3 a /X A + {Q7l) a /X B + {Q7 2 )a/Xc + (Q7 3 )a/x D < 0, (14) 

where xa = A, xb = Vx 2 X 3 > x c = Vx 1 ^ 3 , x ° = ^/x 1 X 2 - The positive parameter 
A has been introduced to put all the variables Q9, Q7i, Q5i, and Q3i on equal footing. 
These supersymmetry conditions can be easily cast in the T-dual form of the Type IIA 
supersymmetry constraints discussed in [2]. 

III. CONSTRUCTIONS OF SM-LIKE STRING VACUA FROM TYPE IIB FLUX 
COMPACTIFICATION 

Similar to the past constructions (see, specifically those of [13]), we construct the SM-like 
models as descendants of the Pati-Salam model based on SU (A) c x SU (2) L x SU (2) R gauge 
symmetry. The hypercharge is 

Qy = Q I3R + ^y I± , (15) 

where the non-anomalous U(1) B -l is obtained from the splitting of the U(A) C branes, i.e. 
U(A) — > U(3)c x U{\)b-l- Similarly, the anomaly-free U(l)i 3R gauge symmetry is from the 
non-Abelian U(2) R or Sp(2) R gauge symmetry. There are three main frameworks to realize 
the Pati-Salam gauge sector in the Type IIB magnetized D-brane scenario (T-dual to the 
Type IIA intersecting D6-brane one): 

(i) The starting observable sector gauge symmetry is U(A) C x U(2) L x U(2) R [13]. In 
this framework, the three "anomalous" gauge symmetries U(\)c, U(1)l and U(1) R can be 
treated as global ones, since the associated gauge bosons obtain masses via B A F Chern- 
Simons couplings. [Those are effective couplings that arise from the D-brane world-volume 
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Chern-Simons couplings and are responsible for the Abelian gauge anomaly cancellation via 
the Green-Schwarz mechanism.] The gauge symmetry breaking chain is of the form: 

SU(A) x SU(2) L x SU(2) R 
-> SU(3) C x SU(2) L x SU(2) R x U(1) B -l 
-> SU(3) C x SU(2) L x C/(1) J3J{ x t/(l) B _ L 

-> 5C/(3) C x 5C/(2) L x C/(l) y , (16) 

where the first and second step can be achieved by splitting the U(A) C branes and U(2) R 
branes in one two-torus direction, and the third step by giving vacuum expectation values 
(VEVs) to the scalar components of right-handed neutrino chiral superfields (or a scalar 
component of an exotic non-chiral chiral superfields) at the TeV scale. Alternatively, we may 
skip the second step and directly break the SU(2) R x U(1) B -l gauge symmetry down to the 
U(l)y by giving VEVs to the scalar components of right-handed neutrino chiral superfields. 
Within this framework one typically obtains enough SM Higgs doublet pairs (from the chiral 
or the non-chiral massless sector) with Yukawa couplings to quarks and leptons and hence 
explain the mass hierarchy and mixings for the fermions in the SM sector at the tree level. 
However, it should also be noted that in a few cases all the SM Higgs doublet pairs have 
the global U(l) quantum numbers that do not allow for Yukawa couplings to quarks and 
leptons; these models therefore face serious phenomenological obstacles. 

(ii) The starting gauge symmetry is one-family U(4) x Sp(8)i x Sp(8) R or two-family 
U (4) x Sp(A) L x Sp(A) R , which can be broken down to the four-family C/(4) x U(2) L x U (2) R 
or £7(4) x SU(2) L x SU(2) R by parallel splitting the D-branes, originally positioned on the 
O-planes, in three or two-tori directions, respectively. [Both the string theory and field 
theory aspects of the brane splittings in this framework are discussed in detail in [14] for 
constructions without fluxes. The flux vacua of that type (with one unit of quantized flux) 
was constructed in [37].] In the field theory picture ("Higgsing") the four-families (/ = 4) 
are obtained when we decompose the original chiral supermultiplets (4,8,1) and (4,1,8), 
or (4,4,1) and (4,1,4) into four copies of (4,2,1) and (4,1,2) after the gauge symmetry 
breaking. These Higgsings, as discussed in [14], preserve the D- and F-flatness, and thus 
the symmetry breaking can take place at the string scale. Thus, in these cases the resulting 
spectrum is that of four-dimensional N — 1 supersymmetric four-family Pati-Salam models. 
The symmetry breaking chains for these two pictures are given respectively by 

SU(A) x Sp(8) L x Sp(8) R 
-> SU(4) x U(2) L x U{2) R 
-> SU(3) C x U(2) L x U(2) R x U{1) B -l 

-> SU(3) C x U(2) L x U{l) Y . (17) 
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and 



SU(4) x Sp(4) L x Sp(4) R 
-> 517(4) x 5t/(2) L x SU(2) R 
-> 5C/(3) C x 5C/(2) L x Sl/(2) fl x 17(1) B -l 

-> 5C/(3) C x 5C/(2) L x U(1) Y ■ (18) 

The first and second step can be achieved by splitting the Sp- and [7(4)-branes at the 
string scale, and the third step by giving VEVs to the scalar components of the right- 
handed neutrino superfields at the TeV scale. Note that for the model with the original 
£7(4) x Sp(8)i x Sp(8) R symmetry, the resulting U(1)l and U{1)r are not anomalous since 
they are part of the non-Abelian Sp symmetries. One expects that at least the gauge boson 
of U(1)l will obtain a mass at the electroweak scale, which is excluded by experiments. In 
order to evade this problem, we allow only for the Sp(8) R gauge symmetry in the SM sector, 
while the Sp(8)l gauge symmetry is not. Moreover, we consider the variants of the U(4) x 
Sp(4)l x Sp(4) R model, i.e. the models with the gauge symmetry 17(A) x U(2)l x Sp(4) R 
or U(4) x Sp(4)l x U(2) r . [This analysis can also be applied to three-family models with 
gauge symmetry U(A) x Sp(6)l x Sp(Q) R where the Sp(6)l x Sp(6) R gauge symmetry can 
be broken down to the Sp(2)l x Sp(2) R by the Higgs mechanism, however this symmetry 
breaking pattern breaks supersymmetry [14], and thus may only be implemented within the 
framework of supersymmetry breaking at scale larger than the electroweak scale.] 

(hi) The starting symmetry is the Pati-Salam-like U(4) c x Sp{2) L x Sp{2) R . Without 
fluxes, the first models of that type were toroidal orientifolds with intersecting D6-branes [18] 
where the RR tadpoles were not explicitly cancelled. The Z 2 x Z 2 orientifold construction in 
Ref. [14] cancelled the RR-tadpoles by introducing an additional stack of branes with unitary 
symmetry in the hidden sector. In the T-dual Type IIB picture those are the magnetized 
D9-branes with a negative D3 charge. In Ref. [36] these types of magnetized D9-branes 
were employed to find the MSSM-like model with one unit of quantized flux turned on. 
In this framework, the starting gauge symmetry is similar to the framework (i); the initial 
framework gauge symmetry U(4)c x Sp(2)l x Sp(2)n can be broken down to the SM gauge 
symmetry as 

517(4) x SU{2) L x SU(2) R 
-> SU(3) C x SU(2) L x SU(2) R x U(1) B _ L 

-> SU(3) C x SU(2) L x U(1) Y , (19) 

where the first step can be achieved again, by splitting the U (4) D-branes at string scale and 
the second step is achieved by giving VEVs to the scalar components of the right-handed 
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neutrino chiral superfields at the TeV scale. Note that in the present case there is no U(1)l- 
However, in the models constructed in Ref. [18, 36], there exists only one SM Higgs doublet 
pair; in this case a generic problem is that only the third family can obtain the tree-level 
masses and it is difficult to give masses to the first two families at the quantum level [18]. 
In addition, the SU(2) R x U(\)b-l gauge symmetry can be broken down to the U(l)y only 
by giving VEVs to the scalar components of the right-handed neutrino chiral superfields at 
the TeV scale in this kind of models. 

The Pati-Salam-type models with only one SM Higgs pair suffer from serious phenomeno- 
logical problems [48]. Even though one may be able to generate the most general Yukawa 
couplings via radiative corrections, the mass matrix of up-type quarks is proportional to 
that of down-type quarks, and the mass matrix of neutrinos is proportional to that of lep- 
tons. (Note the renormalization group equation running effects on these mass matrices are 
negligible.) Therefore, on the one hand, the masses for the quarks, leptons and neutrinos 
satisfy 

m u : m c : m t = m d : m s : m d , m Ve : m Ufi : m Vr = m e : m M : m T . (20) 

The above fermion mass relations are obviously wrong from the known experiments. On 
the other hand, the CKM quark mixing matrix and the PMNS neutrino mixing matrix 
are proportional to the identity matrix which implies that the quark and neutrino mixing 
angles vanish, again in contradiction with experiments. Note that these problems for the 
fermion masses and mixings cannot be solved by loop corrections because the SU(2) R x 
U(1)b-l gauge symmetry is broken at TeV scale. There is a consensus that the minimal 
super symmetric Pati-Salam or Left-Right model should have at least two SM Higgs doublets 
[42]; therefore, the construction (without fluxes) in [14], which actually contains two SM 
Higgs doublets is the first one which can realize the embedding of the supersymmetric 
Pati-Salam model with realistic features in the Type IIA intersecting D-brane scenario, or 
equivalently, the T-dual Type IIB magnetized D-brane scenario. 

The presence of fluxes further complicates the constructions of these types of models: in 
the Type IIB background, the G 3 fluxes give a large positive contribution to the D3-brane 
RR tadpoles, thus making it extremely hard to satisfy the D3 charge tadpole cancellation 
conditions by the magnetized D-brane sectors. In the first model with one unit of quan- 
tized flux [36], the large positive contribution to D3 charges from the flux, is cancelled by 
the "hidden sector" magnetized D9-branes, carrying negative D3 charges (first introduced 
for vacua without fluxes in [14]). Four-family models with one unit of flux and the starting 
SM-sector gauge symmetry (ii) were constructed by introducing a single stack of magnetized 
D9-branes with the negative D3 charge in [37]. These very few specific semi-realistic con- 
structions are extremely constrained, thus implying that semi-realistic flux vacua are hard 
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to come by. In this paper we advance this program in a new direction, by introducing mag- 
netized D9-branes with negative D3 charges as a part of the SM-sector within frameworks (i) 
and (ii). As a consequence we obtain a large number of the three- and four-family SM-like 
flux with as much as three units of flux turned on. In particular these constructions provide 
first four- dimensional N = 1 super symmetric SM-like string vacua (i.e. three units of flux) 
as well as first examples of semi-realistic SM-like string vacua with two units of flux and 
many new models with one unit of flux. In addition, we also obtained SM-like flux models 
where the "hidden sector" D7-brane gauge dynamics can induce gaugino condensation that 
can stabilize the third toroidal Kahler modulus (two other toroidal Kahler moduli are fixed 
by supersymmetry conditions in the D-brane sector); this is an explicit construction that 
may realize the KKLT mechanism [45]. 

In the following subsections we shall present explicit representative models within its 
class. Within each class there are typically more models and a sizable number of models 
within each class has been obtained by running computer code. The representative model 
in each class is typically chosen to have a minimal number of chiral (mainly right) exotics. 
In the following, we give a concise description of the representative models. (Please, refer to 
the Appendix for tables containing these models and a detailed explanation of the notation 
employed.) A concise discussion of phenomeno logical implications of these models will be 
given in [43]. 

A. Models with Supersymmetric Fluxes 

In this subsection, we construct SM-like string vacua with the supersymmetric flux con- 
figuration, i.e. three units of quantized flux. Again, the key feature is the introduction of 
magnetized D9-branes with the negative D3 charge, which is a part of the SM-sector. These 
are the first three- and four-family SM-models within the supersymmetric flux background. 
The D-brane configurations of Model — Ti — 3 and Model — F± — 3 are given in Table VII and 
Table XIV, respectively. The chiral spectrum for the three-family model (Model — T\ — 3) 
is given in Table II. All three toroidal Kahler moduli in these models are fixed by the su- 
persymmetry conditions for the D-brane sector [Note that the open string moduli and the 
Kahler moduli can form combined D-flat directions, corresponding to the brane recombina- 
tion (see Ref. [2] in the context of specific orientifold compactifications). However, due to 
the flux back-reaction it is expected that the open string moduli could become massive [30] ; 
in this case the supersymmetry conditions do stabilize Kahler moduli.] But, the SM Higgs 
doublets do not have Yukawa couplings to quarks and leptons due to the "wrong" quantum 
numbers under the global U(l) symmetries, and one has to look for new ways to generate 
quark and lepton masses. We shall further discuss the masses of the SM families as well as 
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TABLE II: The chiral spectrum in the open string sector of Model — T\ — 3. 
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SM chiral exotics in [43]. 

B. Models with Non-supersymmetric Fluxes 

In this subsection, we shall consider the string vacua with non-supersymmetric fluxes, i.e. 
with two and one units of quantized fluxes. With fewer units of quantized flux, there is more 
freedom in satisfying the tadpole conditions. In the following we only present some typical 
three- and four-family models for each phenomenologically interesting case. Since the gauge 
symmetry breaking chain for each model can be easily determined from the analysis at the 
beginning of this section, we mainly focus on additional phenomenological aspects of these 
models. 

1. Two Flux Units Models 

We constructed the first three- and four-family SM-like string vacua with two units of 
quantized flux with the representative models: Model — Ti — 2 (Table VIII), Model — Fx — 2 
(Table XV), Model - F 2 - 2 (Table XVI) and Model - F 3 - 2 (Table XVII). 

For the three-family model (Model — T x — 2), its chiral spectrum is given in Table III. Note 
that there is one left-chiral exotic (4, 2, 1, 1), which has Yukawa couplings to the right-chiral 
ones and SM Higgs doublet pairs (Higgs bidoublets), thus these exotics can obtain a mass 
at the electroweak scale. In addition, there are five pairs of (non-chiral) SM Higgs doublet 
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pairs, arising in the be sector when the b and c stacks of D-branes are coincident on the 
first two-torus; these SM Higgs doublet pairs have correct global U(l) quantum numbers, 
allowing for the Yukawa couplings to quarks and leptons and thus the SM fermion masses 
and mixings can be generated at the tree level. 



TABLE III: The chiral spectrum in the open string sector of Model — T\ — 2. 
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For the four-family models, there do not exist any left chiral exotics. Similar to the 
Model — 7\ — 2, the suitable SM fermion masses and mixings can be obtained due to the 
Yukawa couplings of the SM Higgs doublet pairs in Model — F 3 — 2. However, the SM 
fermion masses and mixings can not be generated at the tree level in the Model — Ti — 2 and 
Model — F 2 — 2 because the SM Higgs doublet pairs have wrong quantum numbers under 
the U{1) global symmetries and thus no Yukawa couplings to quarks and leptons. 

2. One Flux Unit Models with [7(2)^^ Negative D3 Charge Branes 

When one unit of flux is turned on, there is a wealth of models and these constructions 
of three- or four-family SM-like flux vacua can be classified in the following way: 

(1) Model -Ti-1 (Table IX) and Model -Fx -I (Table XVIII). Except for some symmet- 
ric and anti-symmetric representations, these two models do not contain any bifundamental 
chiral exotics in the observable sector. Their chiral spectra are given in Table IV and Table 
V, respectively. In particular, even though for the four-family model (Model — Ti — 1) there 
is an anomaly-free U(1)r gauge symmetry, it is broken at the "right-handed" scale, when 
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the SU(2) R x U(1)b-l gauge symmetry is broken down to the U(l)y by giving VEVs to 
the scalar components of right-handed neutrino chiral superfields. However, in these two 
models, the SM Higgs doublet pairs do not have Yukawa couplings to quarks and leptons, 
due to the wrong global U(l) quantum numbers. 



TABLE IV: The chiral spectrum in the open string sector of Model — T\ — l. 
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TABLE V: The chiral spectrum in the open string sector of Model — F\ 


- 1. 
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(2) In the Model - T 2 - 1 (Table X), Model - T 3 - 1 (Table XI) and Model - F 3 - 1 
(Table XX), there are four or five pairs of non-chiral Higgs bidoublets arising from the be 
sector which can couple via Yukawa couplings to quarks and leptons, which can generate 
fermion mass hierarchies and mixings. Some of the additional chiral exotics can obtain large 
masses by coupling to these non-chiral Higgs bidoublets. 
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Even though these models employ non-chiral Higgs bidoublets to give masses to fermions, 
we can easily find models with chiral Higgs bidoublets with appropriate Yukawa couplings 
to quarks and leptons, see, e.g., Model — F\ — 2 (Table(XV)). 

Another interesting four-family model is Model — F 4 — 1 (Table XXI); it does not contain 
any chiral exotics in the observable sector. The chiral Higgs bidoublets in the be sector allow 
for the Yukawa couplings to quarks and leptons and thus the fermion masses and mixings 
can be generated at the tree level. Note also that the U(A) C symmetry emerges by Sp(16) 
D7-brane splitting at the string scale and thus U(l)c gauge symmetry is non-anomalous 
and is broken at the "right-handed" scale by the VEVs of the scalar components of the 
right-handed neutrino chiral superfields. 

3. One Flux Unit Models with U(4)c Negative D3 Charge Branes 

In this construction the U(A) C is due to the negative D3 charge magnetized D9-branes. 
Model - T 4 - 1 (Table XII) and Model - F 5 - 1 (Table XXII) are such three- and four-family 
SM-like models. In both models, the SU(2) R gauge symmetry is generated by D7-brane 
splitting, which yields eight and four copies of right-chiral representations, respectively. In 
particular, the four-family model (Model — F 5 — 1), whose chiral spectrum is given in Table 
VI, have several nice phenomenological features: 

(i) No additional bifundamental chiral exotics in the observable sector; 

(ii) The suitable fermion masses and mixings can be generated by tree level Yukawa 
couplings to the SM Higgs doublet pairs; 

(iii) No additional electroweak scale U(l) symmetry. 



TABLE VI: The chiral spectrum in the open string sector of Model — F$ — 1. 
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C. Models with Kahler Moduli Stabilized by D7-brane Gauge Dynamics 

In this subsection, we consider the possibility of realizing a stabilization of the toroidal 
Kahler moduli stabilization a la KKLT mechanism [45]. In the original paper, the KKLT 
vacua are achieved via three steps: 

(1) Turning on self-dual three-form fluxes on Type IIB Calabi-Yau manifold, the flux- 
induced superpotential will fix the dilaton and all the complex structure moduli; 

(2) Introducing a Kahler-moduli dependent non-perturbative superpotential, due to D7- 
brane strong infrared gauge dynamics or Euclidean D3-brane instanton effect. The vacuum 
is a supersymmetric anti-de Sitter one with the Kahler moduli fixed; 

(3) Adding a set of anti-D3-branes to lift the anti-de Sitter vacuum to a de Sitter one. 
We will only focus on the second step, i.e. the generation of the non-perturbative super- 
potential, since it may help us stabilize all the toroidal moduli fields. 

In our framework, the Type IIB SM-like flux vacua typically require a stack of (filler) D3- 
branes, sitting on the 03-planes, in order to cancel the D3 charge tadpoles due to redundant 
negative D3 charges introduced by the magnetized D9-branes. This stack of D3-branes has 
typically a negative beta function; it thus possesses a non-perturbative gauge dynamics that 
results in a non-perturbative superpotential, due to gaugino condensation. However, this 
superpotential depends only on the dilaton-axion field, and is independent Kahler-moduli. 
Thus, in order to generate the non-perturbative superpotential that depends on the toroidal 
Kahler moduli, the strong gauge dynamics has to arise due to D7-branes [45]. 

Recently, it was suggested [39] that the flux-induced soft mass terms may help decouple 
the open-string moduli on D7-branes, leaving an infrared gauge theory with strong dynamics 
on the world- volume of this stack of D7-branes (in the hidden sector). However, in the 
concrete constructions of SM-like flux vacua a stack of hidden sector D7-branes typically do 
not have a negative beta function; this is due to the fact that the magnetized D9-branes 
with negative D3 charge have large "intersecting numbers" with the D7-branes and thus a 
large number of chiral matter, charged under D7-brane gauge symmetry. Additional chiral 
matter typically drastically modifies the infrared gauge dynamics. [Light chiral matter 
can influence the (supersymmetric) gauge dynamics in two key aspects (see, e.g., [46] and 
references therein): (i) the chiral matter contributes to the beta function and thus affect the 
infrared dynamics (phase structure); (ii) it may induce matter condensation and contribute 
to the non-perturbative superpotential.] 

In this paper, we present the first SM-like flux vacua with strong gauge dynamics, re- 
sulting in gaugino condensation, on D7-branes of the hidden sector: Model — T 5 — 1 (Table 
XIII) and Model — F e — 1 (Table XXIII). These models have three- and four-family fermions 
(and additional chiral exotics), respectively. In these models the two (out of three) toroidal 
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Kahler moduli have been fixed by supersymmetry conditions. For the four-family model 
(Model — F 6 — 1), its hidden sector is composed of two stacks of D7-branes denoted by 
(D7)i and (D7) 2 , respectively. Both of them carry Sp(A) gauge symmetry and the associ- 
ated beta functions are -3(0) and -5 (-2), respectively. Here the beta functions in the brackets 
include the one-loop contribution from the open-string moduli. These open-string moduli 
are expected to become massive due to the flux back- reaction (see, e.g., [30]), and in this 
case the strong gauge dynamics can in principle induce the non-perturbative superpotential. 
Note, however, that gauge dynamics of the (_D7)i-branes results in the superconformal or 
Coulomb phase regime and the non-perturbative superpotential cannot be dynamically gen- 
erated. As a result, only (D7)2-branes can generate a Veneziano-Yankielowicz-type superpo- 
tential, induced by gaugino and matter condensations. A similar analysis can be applied to 
the three-family model (Model — T 5 — 1): in the case that open string moduli are decoupled 
(the beta function is -9(0)), the Veneziano-Yankielowicz superpotential is generated due to 
the gaugino and matter condensates on the D7-brane. 

IV. DISCUSSIONS AND CONCLUSIONS 

In this paper we have advanced a program for explicit constructions of supersymmetric 
SM-like string vacua with supergravity fluxes turned on. In particular we obtained large 
classes of such Type IIB models on Z 2 x Z 2 orientifolds with one, two, and three units of 
quantized flux turned on. These models provide an important stepping stone toward broader 
classes of realistic constructions that not only contain the three- (or four-) family SM sector, 
but also stabilize some of the moduli (typically all toroidal moduli can be fixed). 

Before our work, techniques for constructions of the chiral D-brane sectors flux vacua on 
Type IIB orientifolds were developed in [30, 31], however, the semi-realistic constructions 
remained elusive until recently. The technical reason for such difficulties are large positive 
three-form flux contributions to the D3 charge in the internal space (D3-branes RR-tadpole), 
which makes the D3 charge conservation constraint hard to satisfy. In the first examples 
of SM flux vacua on Z 2 x Z 2 orientifolds this constraint has been satisfied [36, 37] by the 
introduction of the hidden sector magnetized D9-branes which carry negative D3 charges. 
[This type of D-branes were first introduced in the Type IIA context with intersecting 
D6-branes without fluxes in [14].] Within this framework one three-family [36] and one 
four-family [37] SM-like models with one unit of quantized flux were obtained. However, 
one remains to be confronted by serious problems: 

(i) There remains the outstanding problem of constructing semi-realistic SM-like string 
vacua with supersymmetric fluxes, that for Z 2 x Z 2 orientifolds correspond to three units 
of quantized flux. The focus therefore shifted to the study of non-supersymmetric flux 
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effects as the key to break supersymmetry, and a detailed study of soft supersymmetry 
breaking masses due to fluxes [39, 40]. The flux-induced supersymmetry breaking, however, 
leads typically to soft masses M so f t ~ j^-, which implies an intermediate string scale or 
inhomogeneous warp factor in the internal space in order to stabilize the electroweak scale; 

(ii) In spite of the successes in stabilizing the dilaton and toroidal complex structure 
moduli, the toroidal Kahler moduli are not completely fixed, thus one remains faced with 
the large vacuum degeneracy problem; 

(iii) The fact that (supersymmetric) supergravity fluxes are abundant, makes it is im- 
perative to extend the constructions to semi-realistic models with more than one unit of 
quantized flux. 

In this paper, we have made important progress in addressing the above issues. The 
key ingredient in the new SM-like flux constructions is the introduction of the negative D3 
charge magnetized D9-brane as a part of the SM sector. These constructions turned out 
to be less constraining and resulted in three- and four-family SM-like string vacua with up 
to three units of quantized flux, thus leading to the first fully supersymmetric SM-like flux 
vacua. In addition to toroidal complex structure moduli and the dilaton being fixed by the 
flux, the supersymmetry conditions in the D-brane sector typically fix all the toroidal Kahler 
moduli, and the string scale is close to the Planck scale. We also constructed the first three- 
and four-family SM-like models with two units of the quantized flux and many new models 
with one unit of the quantized flux. Typically the representative models have (mainly right-) 
chiral exotics in the SM sector; however, we have also presented a few three- and four-family 
models with no SM-sector chiral exotics. (Note, the models do have additional tensor fields, 
i.e. chiral superfields in the symmetric and/or anti-symmetric representation of the unitary 
gauge symmetry, typically associated with the negative D3 charge magnetized D9-brane.) 
We have also been able to construct SM-like flux vacua with strong infrared gauge dynamics 
on the hidden sector D7-branes, which leads to non-perturbative superpotential that can fix 
a remaining toroidal Kahler modulus. 

The SM Higgs doublet pairs appear at the intersections of the SU(2) L and SU(2) R . For 
most models, the SM fermion masses and mixings can be generated due to the tree-level 
Yukawa couplings of such SM Higgs doublet pairs to quarks and leptons. However, in a 
few cases, most notably for the three-family SM-like model with the supersymmetric flux, 
all such Higgs fields have the wrong global U(l) quantum numbers that prevent them from 
coupling to quarks and leptons via Yukawa couplings; these models therefore face serious 
phenomenological difficulties. Note also that typically in these models some of the SM chiral 
exotics can obtain masses at least at the electroweak scale, due to the Yukawa couplings of 
the SM Higgs doublet pairs to such exotics. 

In spite of a number of advances made with these new constructions, there are open 
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problems which deserve further study. In particular, further discussion of the specific Yukawa 
couplings to quarks, leptons and chiral exotics, as well as the subsequent further implications 
for the masses and mixings of the SM chiral fermions are needed [43]. For SM-like models 
with supersymmetric fluxes, one has to address the supersymmetry breaking mechanism. 
In principle the supersymmetry breaking could be due to the hidden D-brane strong gauge 
dynamics; however, the present models do not possess such a sector. Complete stabilization 
of all the moduli (and not only the toroidal closed sector one) remains an open problem. 
We postpone these issues for future research. 

Acknowledgments 

We would like to thank Paul Langacker for many discussions, encouragement, and collab- 
oration on topics related to this paper. The research was supported in part by the National 
Science Foundation under Grant No. INT02-03585 (MC) and PHY-0070928 (T. Li), by the 
Department of Energy Grant DOE-EY-76-02-3071 (MC, T. Liu), and by the Fay R. and 
Eugene L. Langberg Chair (MC). 



[1] M. Cvetic, G. Shiu and A. M. Uranga, Phys. Rev. Lett. 87, 201801 (2001), hep-th/0107143. 
[2] M. Cvetic, G. Shiu and A. M. Uranga, Nucl. Phys. B615, 3 (2001), hep-th/0107166. 
[3] R. Blumenhagen, L. Gorlich, B. K6rs and D. Lust, JHEP 0010 (2000) 006, hep-th/0007024. 
[4] G. Aldazabal, S. Franco, L. E. Ibanez, R. Rabadan and A. M. Uranga, JHEP 0102, 047 
(2001), hep-ph/0011132. 

[5] G. Aldazabal, S. Franco, L. E. Ibanez, R. Rabadan and A. M. Uranga, J. Math. Phys. 42, 

3103 (2001), hep-th/0011073. 
[6] R. Blumenhagen, B. K6rs and D. Lust, JHEP 0102 (2001) 030, hep-th/0012156. 
[7] C. Angelantonj, I. Antoniadis, E. Dudas and A. Sagnotti, Phys. Lett. B 489, 223 (2000), 

hep-th/0007090. 
[8] C. Bachas, hep-th/9503030. 

[9] M. Berkooz, M.R. Douglas and R.G Leigh, Nucl. Phys. B 480, 265 (1996), hep-th/960613. 
[10] L. E. Ibanez, F. Marchesano and R. Rabadan, JHEP 0111, 002 (2001), hep-th/0105155; 
R. Blumenhagen, B. K6rs and D. Lust, T. Ott, Nucl. Phys. B616 (2001) 3, hep-th/0107138; 
D. Cremades, L. E. Ibanez and F. Marchesano, Nucl. Phys. B643, 93 (2002), hep-th/0205074; 
D. Cremades, L. E. Ibanez and F. Marchesano, JHEP 0207, 022 (2002), hep-th/0203160; 
D. Bailin, G. V. Kraniotis, and A. Love, Phys. Lett. B 530, 202 (2002); Phys. Lett. B 547, 43 



22 



(2002); Phys. Lett. B 553, 79 (2003); JHEP 0302, 052 (2003); C. Kokorelis, JHEP 0209, 029 
(2002); JHEP 0208, 036 (2002); hep-th/0207234; JHEP 0211, 027 (2002); hep-th/0210200; 
hep-th/0406258; hep-th/0412035; T. Li and T. Liu, Phys. Lett. B 573, 193 (2003), hep- 
th/0304258; 

[11] M. Cvetic and I. Papadimitriou, Phys. Rev. D 68, 046001 (2003), hep-th/0303083; Phys. Rev. 

D67, 126006 (2003), hep-th/0303197; M. Cvetic, I. Papadimitriou and G. Shiu, Nucl. Phys. 

B659, 193 (2003) hep-th/0212177; 
[12] R. Blumenhagen, L. Gorlich and T. Ott, JHEP 0301, 021 (2003), hep-th/0211059; G. Ho- 

necker, Nucl. Phys. B666, 175 (2003), hep-th/0303015; G. Honecker and T. Ott, hep- 

th/0404055. 

[13] M. Cvetic, T. Li and T. Liu, Nucl. Phys. B 698, 163 (2004), hep-th/0403061. 
[14] M. Cvetic, P. Langacker, T. Li and T. Liu, hep-th/0407178. 

[15] M. Cvetic, P. Langacker and G. Shiu, Phys. Rev. D 66, 066004 (2002), hep-ph/0205252; 

M. Cvetic, P. Langacker and G. Shiu, Nucl. Phys. B 642, 139 (2002), hep-th/0206115; 
[16] R. Blumenhagen, D. Lust and S. Stieberger, JHEP 0307, 036 (2003), hep-th/0305146. 
[17] D. Lust and S. Stieberger, hep-th/0302221. 

[18] D. Cremades, L.E. Ibahez and F. Marchesano, JHEP 0307, 038 (2003), hep-th/0302105. 
[19] M. Cvetic and I. Papadimitriou, Phys. Rev. D 68, 046001 (2003) [Erratum-ibid. D 70, 029903 
(2004)], hep-th/0303083. 

[20] S. A. Abel and A. W. Owen,Nucl. Phys. B 663, 197 (2003), hep-th/0303124; Nucl. Phys. B 

682, 183 (2004), hep-th/0310257; S. A. Abel and B. W. Schofield, hep-th/0412206. 
[21] D. Lust, P. Mayr, R. Richter and S. Stieberger, Nucl. Phys. B 696, 205 (2004), hep-th/0404134; 

G. L. Kane, P. Kumar, J. D. Lykken and T. T. Wang, hep-ph/0411125. 
[22] R. Blumenhagen, JHEP 0311, 055 (2003), hep-th/0310244; I. Brunner, K. Hori, K. Hosomichi 

and J. Walcher, hep-th/0401137; R. Blumenhagen and T. Weigand, JHEP 0402, 041 (2004), 

hep-th/0401148; 

[23] T. P. T. Dijkstra, L. R. Huiszoon and A. N. Schellekens, hep-th/0403196; T. P. T. Dijkstra, 

L. R. Huiszoon and A. N. Schellekens, hep-th/0411129. 
[24] M. Cvetic, P. Langacker and J. Wang, Phys. Rev. D 68, 046002 (2003), hep-th/0303208. 
[25] S. Gukov, C. Vafa and E. Witten, Nucl. Phys. B 584, 69 (2000) [Erratum-ibid. B 608, 477 

(2001)], hep-th/9906070. 
[26] K. Behrndt and M. Cvetic, hep-th/0403049. 
[27] K. Behrndt and M. Cvetic, hep-th/0407263. 
[28] D. Lust and D. Tsimpis, hep-th/0412250. 

[29] B. S. Acharya, F. Denef, C. Hofman and N. Lambert, hep-th/0308046. 



23 



[30] J.F.G. Cascales and A.M. Uranga, JHEP 0305, Oil (2003), hep-th/0303024. 
[31] R. Blumenhagen, D. Lust and T.R. Taylor, Nucl. Phys. B 663, 319 (2003), hep-th/0303016. 
[32] S. B. Giddings, S. Kachru and J. Polchinski, Phys. Rev. D 66, 106006 (2002), hep-th/0105097. 
[33] S. Kachru, M. B. Schulz, P. K. Tripathy and S. P. Trivedi, JHEP 0303, 061 (2003), hep- 
th/0211182. 

[34] A. Giryavets, S. Kachru and P. K. Tripathy, JHEP 0408, 002 (2004), hep-th/0404243; O. De- 
Wolfe, A. Giryavets, S. Kachru and W. Taylor, hep-th/0411061. 
[35] A. Font, JHEP 0411, 077 (2004), hep-th/0410206. 
[36] F. Marchesano and G. Shiu, hep-th/0408059; hep-th/0409132. 
[37] M. Cvetic and T. Liu, hep-th/0409032. 
[38] C. Kokorelis, hep-th/0309070. 

[39] P. G. Camara, L. E. Ibahez and A. M. Uranga, Nucl. Phys. B 689, 195 (2004, hep-th/03 11241; 

L. E. Ibahez , hep-ph/0408064; P. G. Camara, L. E. Ibahez and A. M. Uranga, hep-th/0408036; 

A. Font and L. E. Ibahez , hep-th/0412150. 
[40] D. Lust, S. Reffert and S. Stieberger, hep-th/0410074; D. Lust, S. Reffert and S. Stieberger, 

hep-th/0406092. 
[41] F. Marchesano, G. Shiu and L. T. Wang, hep-th/0411080. 

[42] See, e.g., "Unification and Supersymmetry: The Frontiers of Quark-Lepton Physics" by 

R. N. Mohapatra; New York : Springer- Verlag, c2003. 
[43] M. Cvetic, P. Langacker, T. Li and T. Liu, UPR-1107-T, in preparation. 
[44] S. Kachru, M. B. Schulz and S. Trivedi, JHEP 0310, 007 (2003), hep-th/0201028. 
[45] S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, Phys. Rev. D 68, 046005 (2003), hep- 

th/0301240. 

[46] V. Kaplunovsky and J. Louis, Nucl. Phys. B 422, 57 (1994), hep-th/9402005; K. A. Intriligator 

and P. Pouliot, Phys. Lett. B 353, 471 (1995), hep-th/9505006. 
[47] E. Witten, Phys. Lett. B 117, 324 (1982). 

[48] We thank Paul Langacker for extensive discussions on phenomenological implications of these 
types of models. See also [43] . 



24 



Appendix: D-brane Configurations and Intersection 
Numbers for SM-like Flux Vacua 

In this Appendix, we tabulate D-brane configurations and intersection numbers for the 
representative three- and four-family models within our new setup. In the first column of 
each table, a, b and c denote the £7(4) (£p(16)), U{2) L (Sp(2) L or £p(4) L ), and U(2) R 
(Sp(A) R or Sp(8)u) stacks of branes, respectively. D3, (D7)i, (D7) 2 , and (D7) 3 represent 
the filler branes along respective QR, QRw, VtROuo and QR9 orientifold planes, resulting in 
Sp(N) gauge groups. N, in the second column, corresponds to the number of filler D-branes 
in each stack. The third column depicts the wrapping numbers of the various D-branes. 
The intersection numbers between the various D-brane stacks are given in the remaining 
right columns where b' and d are respectively the QR images of b and c.In addition, the 
number of symmetric and antisymmetric chiral superfield representations for specific D- 
brane configurations is given. For convenience, we also tabulate the relations among the 
three toroidal Kahler moduli parameters Xii imposed by the supersymmetry conditions. 
The model labels " Model — (T, F)i — n" appearing in the tables denote the " Model — 
(three, four family)i — (n units of fluxes) ". [Since all the models have an even number 
of chiral supermultiplets in the fundamental representation of the Sp(N) gauge groups, these 
models are automatically free of discrete global gauge anomalies [47] .] Finally, we emphasize 
that in this paper we do not fix the convention between the chirality and the sign of the 
intersection number. Instead, we consider the £77(2) D-branes that carry the more realistic 
chiral spectrum (typically only three- or four- families) as the SU(2)l D-branes. These 
representative models therefore possess (mainly) right- chiral exotics. 



TABLE VII: D-brane configurations and intersection numbers for Model — T\ — 3. 



Model - Ti 


- 3 




[U{A) C x U(2) L x 


U (2) n]observable 




j 


N 


(n , m 1 )(n 2 , m 2 )(n 3 , to 3 ) 


"m 




b 


V 


c 


d 


Kahler moduli 


a 




8 


(1,0)(1,1)(1,-1) 








-3 


1 


12 


-10 


X3 = X2 = 2xi 


b 




4 


(1,1)(2,-1)(1,0) 


-2 


2 






6 


-6 


Xs^vTO 


c 




4 


(-2,-l)(4,l)(3,l) 


-46 


-146 













25 



TABLE VIII: D-brane configurations and intersection numbers for Model — T\ 



Model 


-Ti- 


2 


[C/(4) C X U(2) L X J7(2) ii ] 0hs ert,aM e X [Sp(A)] H idden 


j 


N 


(n 1 , m 1 ) (n 2 , m 2 ) (n 3 , m 3 ) 




n R 


b 


V 


c 


d 


a 


8 




(1,0)(1,1)(1,-1) 








-3 


1 


8 


-8 


b 


4 




(2,1)(2,-1)(1,0) 















-4 


c 


4 




(-2,-l)(3,l)(3,l) 


-32 


-112 










(D7) 2 


4 


(0,1)(1,0)(0,-1) 


X3 


= X2 




Xi 




V2T 



TABLE IX: D-brane configurations and intersection numbers for Model — T\ 



Model - Ti - 1 


[C/(4) C X Sp(2) L X U(2) R ]observable X [Sf>( 4 )] Hidden 


j 


N 


(n , m 1 )(n 2 , m 2 )(n 3 , m ) 


"m 


"H 


b 


c 


c' 


a 
b 

c 


8 
2 
4 


(l,0)(3,l)(3,-l/2) 

(0,1)(0,-1)(2,0) 
(-2,-l)(4,l)(3,l/2) 





-23 





-73 


-3 


3 
8 





D'i 


4 


(1,0)(1,0)(2,0) 


-.^ — ^ 12 i 14 -i 

X2 - X3, ^ + X1X2 " 1 



TABLE X: D-brane configurations and intersection numbers for Model — T2 



Model 


-T2-I 


[U{A) C x (7(2) L x [/(2)ij] 


Obseri 


abie 


X 


[Sp(4 


) X Sp(2)] H idden 


j 


N 


(n 1 , m 1 ) (n 2 , m 2 ) (n 3 , to 3 ) 


"m 


"R 


6 


V 


c 


c' 


a 


8 


(l,0)(l,l)(l,-l/2) 








3 


-2 


-4 


4 


b 


4 


(2,-l)(l,0)(5,l/2) 


3 


-3 









-4 


c 


4 


(-2,l)(3,-l)(3,-l/2) 


16 


56 










D3 


4 


(1,0)(1,0)(2,0) 




X3 = 




X2 




|xi = V39 


(D7) 2 


2 


(0,1)(0,-1)(2,0) 
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TABLE XI: D-brane configurations and intersection numbers for Model — T 3 — 1. 



Model 


-T3-I 


[U{A) C X U(2) L X U(2) R } bservable 


X 


[Sp(8)] Hidden 


j 


TV 


(n 1 , m 1 ) (n 2 , m 2 ) (n 3 , to 3 ) 




n R 


6 


b' 


C 


c' 


a 


8 


(1,0)(1,1)(1,-1) 








3 


1 


-3 


3 


b 


4 


(2,-l)(l,0)(l,2) 


-6 


6 









12 


c 


4 


(-2,1)(2,-1)(2,-1) 


10 


54 










(D7) 3 


8 


(0,1)(0,-1)(1,0) 


X3 


= X2 = 


= 3X1 





TABLE XII: D-brane configurations and intersection numbers for Model — T4 — 1. 



Mode; - T 4 - 1 


[J7(4)c X C/(2) L X Sp(4) R } bservable X [Sp(4)] fficWe „ 


j 


TV 


(n , m 1 )(n 2 , m 2 )(n 3 , to 3 ) 


"m 


n R 


6 


7/ 


c 


a 
6 

c 


8 
4 
4 


(-1,-1)(2,1)(2,1) 
(1,0)(3,1)(1,-1) 
(1,0)(0,1)(0,-1) 


-2 
4 



-30 
-4 



3 


-5 


-4 



733 


4 


(1,0)(1,0)(1,0) 


3X3 


= X; 




xi 


+ 8 =1 

X1X2 



TABLE XIII: D-brane configurations and intersection numbers for Model — T5 — 1, here (5 9 are 
beta functions for the associated gauge symmetries in the hidden sector. 



Model 


-T5-I 


[U(4) C X [7(2) L X Sp(8) R ]observable 


X 


[5p(8) X S P (8)} H idden 


j 


N 


(n 1 , to 1 )^ 2 , m 2 )(n 3 , to 3 ) 






b 


b' 


c 


a 


8 


(1,0)(1,1)(1,-1) 








3 


-3 


-1 


b 


4 


(-2,-l)(2,l)(2,l) 


-10 


-54 






-4 


c 


8 


(0,1)(0,-1)(1,0) 














D3 


8 


(1,0)(1,0)(1,0) 




X2 = 


4 _i_ 8 1 

X3 ' 5g + X!X2 1 


(D7) 2 


8 


(0,1)(1,0)(0,-1) 


^ 3 


= -14(-5), /?f D7)a = -9(0) 
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TABLE XIV: D-brane configurations and intersection numbers for Model — F\ 



Model - F x 


-3 




[U(4) C X Sp{A) L X U(2) R } bservab 


e 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 






b 


c 


d 


Kahlcr moduli 


a 
b 

c 


8 
4 
4 


(1,0)(2,1)(1,-1) 
(0,1)(0,-1)(1,0) 
(-2,-l)(4,l)(3,l) 


2 


-46 


-2 

-146 


-2 


8 
8 


-12 


X2 = 2X3 
24 + 20 =1 

X2 X1X2 



TABLE XV: D-brane configurations and intersection numbers for Model — F\ 



Model 


-Fi-2 


[C/(4) C X f/(2) L X U{2) R ] bservable X [S'p(4)] ffidde „ 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 




"H 


& 


V 


c 


d 


a 


8 


(1,0)(2,1)(1,-1) 


2 


-2 


-4 





4 


-10 


b 


4 


(1,1)(2,-1)(1,0) 


-2 


2 






■5 


-3 


c 


4 


(-2,-l)(3,l)(3,l) 


-32 


-112 










(D7) 2 


4 


(0,1)(1,0)(0,-1) 


X2 = 


= 2X3 




2^ 


l 


= 3V6 



TABLE XVI: D-brane configurations and intersection numbers for Model — F2 



Model -F2-2 


[U{A) C X 5 , p(4) L X U(2) R ]observable X [5p(8) X Sp{4)] mdden 


j 


N 


(n 1 , m 1 ) (n 2 , m 2 ) (n 3 , m 3 ) 


"m 


"R 


& 


c 


c' 


a 
b 

c 


8 
4 
4 


(1,0)(2,1)(1,-1) 
(0,1)(0,-1)(1,0) 
(-2,-l)(3,l)(3,l) 


2 


-32 


-2 


-112 


-2 


4 
6 


-10 


D3 
(D7h 


8 
4 


(1,0)(1,0)(1,0) 
(0,1)(1,0)(0,-1) 


2X3 = X2 

4+ 18 -1 

X2 X1X2 
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TABLE XVII: D-brane configurations and intersection numbers for Model — F 3 — 2. 



Model 


-F3-2 


[U{A)c x U(2) L x J7(2)i ? ] 06sert , Q bi e x [Sp{A)\ Hidden 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 






b 


?/ 


c 


c' 


a 


8 


(1,0)(2,1)(1,-1) 


2 


-2 


-3 


-1 


4 


-10 


b 


4 


(2,1)(1,-1)(1,0) 


2 


-2 









-8 


c 


4 


(-2,-l)(3,l)(3,l) 


-32 


-112 










(D7h 


4 


(0,1)(1,0)(0,-1) 


X2 = 


= 2 X3 




5x1 = 


= 3V3 



TABLE XVIII: D-brane configurations and intersection numbers for Model — F\ — 1. 



Model 


-Fi-1 


[U(4) C X U(2) L X Sp(8) R }ob S ervable X [Sp(8)]mdden 


j 


TV 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 


"m 


n R 


b 


?/ 


c 


a 


8 


(1,0)(1,1)(2,-1) 


-2 


2 


4 





-1 


b 


4 


(-2,-l)(2,l)(2,l) 


-10 


-54 






-4 


c 


8 


(0,1)(0,-1)(1,0) 














(D7) 2 


8 


(0,1)(1,0)(0,-1) 


X3 = 


2X2, 


2 

- 


+ 6 =1 

X1X2 



TABLE XIX: D-brane configurations and intersection numbers for Model — F2 — 1. 



Model 


-F2-I 


[t/(4) C X f7(2) L X U{2) R ]observable 


X 


[Sp(W) x 5p(4)] ffid(ie „ 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 




n R 


b 


?/ 


c 


c' 


a 


8 


(1,0)(1,1)(2,-1) 


-2 


2 


4 





-4 


10 


6 


4 


(3,-l)(l,0)(2,l) 


-2 


2 






5 


5 


c 


4 


(-2,1)(3,-1)(3,-1) 


32 


112 










£>3 


16 


(1,0)(1,0)(1,0) 




X3 = 


= 2 X 2 = 


§Xi - \/30 


(£7) 3 


4 


(0,1)(0,-1)(1,0) 
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TABLE XX: D-brane configurations and intersection numbers for Model — F 3 



Model 


-F3-I 


[U{A)c X U(2) L X U(2) R ]observable 




Sp{8)] Hidden 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 






6 


b' 


c 


c' 


a 


8 


(1,0)(1,1)(1,-1) 








-4 


2 


4 


-6 


b 


4 


(2,1)(3,-1)(1,0) 


-2 


2 









4 


c 


4 


(-2,-l)(2,l)(3,l) 


-18 


-78 










(D7h 


8 


(0,1)(1,0)(0,-1) 


X3 = 


= X2 




Ixi 


= V2T 



TABLE XXI: D-brane configurations and intersection numbers for Model — F4 



Model - F4-I 


[5p(16) c x U(2) L x U(2) R ] observaUe 


j 


N 


(n 1 , m 1 )^ 2 , m 2 )(n 3 , m 3 ) 




"H 


6 


b' 


c 


d 


Kahler moduli 


a 
b 

c 


16 

4 
4 


(1,0)(1,0)(1,0) 
(2,-l)(0,l)(3,-l) 
(-2,-l)(4,l)(l,l) 




-10 
-6 




10 

-58 


1 




-1 

64 





X1X3 = 6 
Xl + 12 = 1 

X2 X1X2 



TABLE XXII: D-brane configurations and intersection numbers for Model — F5 



Model 


-F 5 -l 


[U(A) C X J7(2) L X Sp(8) R ]observable X 


[S , p(4)]_ fJ - id(ie „ 


j 


N 


(n 1 , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 




n R 


6 


?/ 


c 


a 


8 


(-1,-1)(2,1)(2,1) 


-2 


-30 


-4 





1 


b 


4 


(1,0)(1,1)(2,-1) 


-2 


2 









c 


8 


(1,0)(1,0)(1,0) 














{Dl)i 


4 


(1,0)(0,1)(0,-1) 


X3 = 


2X2, 3I 


f 3 =1 

X1X2 



30 



TABLE XXIII: D-brane configurations and intersection numbers for Model — Fq — 1, here (3 9 are 
beta functions for the associated gauge symmetries in the hidden sector. 



Model - F 6 - 1 


[U(4) C X Sp(8) L X U(2) R ] 0b servable X [S P (4) X Sp(4)] ffi(Me „ 


j 


N 


(n , m 1 )(n 2 , m 2 )(n 3 , m 3 ) 


"m 




6 


c 


c' 


a 
b 

c 


8 
8 
4 


(1,0)(1,1)(1,-1) 
(0,1)(0,-1)(1,0) 
(-1,-1)(3,1)(2,1) 




-4 





-44 


-1 


G 
3 


-4 


(D7)i 
(D7h 


4 
4 


(1,0)(0,1)(0,-1) 
(0,1)(1,0)(0,-1) 


X2 X3, ^| + X1X3 1 

= -3(0), P 9 {D7)2 = -5(-2) 
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